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I. INTRODUCTION 



This thesis investigated the performance of a low cost, commercial, space config- 
ured, Global Positioning System (GPS) receiver 1 flown aboard Space Shuttle Discovery, 
STS-51, as part of DTO 0700-6, 2 GPS On-orbit Demonstration (GOOD). The DTO was 
sponsored and funded by NASA Johnson Space Center, and developed with the support 
of NASA contractors, and students, including the authors, from the Naval Postgraduate 
School. Data recorded during the mission was analyzed to evaluate navigational perfor- 
mance of the receiver. Additionally, post-flight filtering of this data was performed in 
order to determine whether a significant increase in performance could be obtained 
through filtering. 

A. BACKGROUND 

1. Space Shuttle Orbiter Baseline Navigation System 

A wide variety of equipment is employed in the Orbiter’s baseline navigation 
system. All navigation sensor information is supplied to a six-state suboptimal Kalman 
filter, which provides the navigation functions with three position and three velocity 
states. The three position states are the coordinates specifying the Orbiter’s position vec- 
tor in the Aries-mean-of-1950 (M50) Cartesian coordinate system. 3 Likewise, the three 
velocity states define the Orbiter’s velocity vector in the M50 system. 



Trimble Advanced Navigation Sensor (TANS) Quadrex GPS Receiver Processor Unit. 
"Detailed Test Objective Number 0700-6. 

T 

The M50 system is defined in NASA Technical Memorandum X-58153, October 1974. 
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During the ascent phase of a mission, the Inertial Measurement Unit (IMU) 
is the primary sensor, providing attitude and acceleration data to the Kalman filter. This 
data is augmented by ground based C-band radar-tracking information uplinked over an 
S-band communication link. During the on-orbit coasting phase of a mission, the IMU 
provides attitude data, and acceleration data from Orbital Maneuvering System (OMS) 
burns. Accelerations falling below the IMU threshold arise from opposing Reaction 
Control System (RCS) thrusters that do not form a perfect couple (vernier effect), and 
from external venting of gasses and waste products. These unaccounted for accelerations 
result in steadily increasing navigational error. 

While on-orbit, the ground continues to track the Orbiter using ground-based 
C-band radar. Two-way Doppler tone ranging over S-band and Ku-band communication 
links, either direct or via the Tracking and Data Relay Satellite (TDRS) system, provides 
additional tracking capability. When the Orbiter’s navigational state is observed to devi- 
ate from the ground based tracking trajectory by a pre-defined, mission-dependent 
amount, a new state vector consisting of three position states, three velocity states, and a 
time tag is uplinked by Mission Control. At various times during a mission (prior to ren- 
dezvous and deorbit burn), an IMU alignment may be performed using an on board star- 
tracker to correct attitude error caused by gyro drift. 

During the re-entry phase of a mission, IMU data is augmented with drag 
modeling data from 250,000 feet down. As the Orbiter passes through the ionosphere, 
all radio-navigation and communication signals are blacked out for a period of time. 
Upon exit from blackout, contact is first established by C-band tracking radar. A state 
vector can be uplinked as soon as S-band communication is regained. Subsequently, the 
Orbiter can receive L-band Tactical Air Navigation (TACAN) station signals and begin 
area navigation, combining TACAN range and bearing data with barometric (30,000- 
2,500 feet) and radar (2,500 feet down) altimeter measurements. Final approach and 
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landing are accomplished with a microwave scan beam landing system (MSBLS), begin- 
ning normally at 10,000 feet, 10 nautical miles downrange from touchdown. 

2. Space Shuttle Orbiter GPS Navigation System 

Pursuant to a study contract commissioned by NASA, Rockwell Internation- 
al’s Space Systems group conducted a design and integration study of a GPS-based pri- 
mary navigation system for the Orbiter in the late 1970’s. (Van Leeuwen et al, 1979, pp. 
118-135) The study demonstrated that the use of on board satellite GPS receivers for 
precise orbit determination was clearly feasible, and expected the improved navigation 
capabilities to yield significant operational benefits. The study concluded the GPS-based 
system to be a technically sound and cost-effective proposition. Based on the study, 
plans were laid to install a GPS navigation system in all Shuttle orbiters beginning in the 
early 1980’s, with follow-on goals of deleting certain equipment from the baseline navi- 
gation system. Within about two years, however, the decision to install GPS was 
reversed in favor of continuing with the baseline navigation configuration. Certain GPS 
provisions, notably antennas, cabling, and bulkhead feedthroughs, were nevertheless 
retained, and currently exist on all Orbiter vehicles. (Saunders, 1994, pp. 1-13) 

The issue of GPS installation in the Orbiter fleet surfaced again in the early 
1990’s. Renewed interest was motivated by the planned phase out of TACAN stations. 
Since TACAN was used as the Orbiter’s primary navigation aid following exit from 
blackout, through MSBLS acquisition, NASA considered suitable alternatives. Looking 
at the direction the Department of Defense (DOD) and the Federal Aviation Administra- 
tion were heading in, GPS was chosen as the replacement for TACAN. A developmental 
test for the Orbiter GPS navigation system flew aboard STS-61 in December 1993, and 
the system is presently expected to be operational in 1996. (Kachmar et. al., 1993, pp. 
313-326) 
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3. GPS On-orbit Demonstration 

In mid- 1992, with the foundation for installing an Orbiter GPS navigation 
system laid, the crew of STS-51 conceived the GOOD DTO as a low cost pathfinder 
project, to look at GPS in orbit, to quantify advantages, and identify potential problem 
areas for Space Shuttle operations. Data from the DTO could then be used to comple- 
ment the more highly integrated GPS Development Flight Test. Since STS-51 would 
carry another payload with its own GPS receiver, the Orbiting Retrievable Far and 
Extreme Ultraviolet Spectrometer-Shuttle Palette Satellite (ORFEUS-SPAS), the DTO 
would also permit the evaluation of relative GPS 1 . Successful utilization of GPS on the 
Orbiter could show benefits for use on other programs, such as Space Station, or for use 
as a utility with other primary and secondary payloads, which require precise location 
and timing information. Initially, goals of this experiment were as follows: 

• Evaluate receiver performance in orbit by comparing its state vector to that deter- 
mined by ground tracking and Orbiter IMU’s. 

• Evaluate the number and location of GPS antennas required to provide best naviga- 
tion solutions for flight deck experiment applications. 

• Determine the quality of GPS data received during on-orbit operations by collecting 
GPS health data. 

• Evaluate the accuracy of relative GPS, using GPS receivers both in the crew cabin 
and on ORFEUS-SPAS, with Orbiter radar and laser rangefinders as a reference. 

• Evaluate postflight the accuracy of relative GPS using data from Orbiter and SPAS 
GPS receivers. 

One of the computer displays developed for this DTO showed the magnitude 
of the position difference between the Orbiter GPS and Orbiter IMU based state vectors 
versus time. 



1 Aspects of relative GPS are covered in the thesis “Theoretical Basis for State Vector Compari- 
son, Relative Position Display, and Relative Position/Rendezvous Prediction” by LT Lester Makepeace, 
and the thesis “NPS State Vector Analysis and Relative Motion Plotting Software for STS-51” by LT Lee 
Barker. 
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Since error in the IMU based state vector increased with time, the root sum 
square (RSS) difference (delta) between the GPS state vector and IMU state vector was 
expected to increase with time. This difference was expected to collapse to zero when 
Mission Control uplinked a new state vector based on ground tracking. This behavior 
was first observed in orbit on flight day three, when a 20,000 ft. delta collapsed to about 
300 feet following an update. An illustration of a similar event on the computer display 
is shown in Figure 1 (the x-axis represents time, the y-axis represents RSS difference). 
The close correlation between expected behavior and actual behavior indicated the GPS 
position solution to be near truth. 




Figure 1: State Vector Differences vs. Time 
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B. GPS OVERVIEW 



The GPS constellation consists of 21 operational satellites, and three active spares, 
distributed in six orbital planes with three or four operational satellites in each plane. 
The ascending nodes of each plane are separated by 60° intervals, and each plane has an 
inclination relative to the equator of 55°. The satellites orbit at an altitude of 20,200 km, 
with a corresponding period of 12 hours. In comparison, the Space Shuttle orbits at an 
altitude of approximately 300 km, with a corresponding period of l'A hours. The satel- 
lites are positioned so that a minimum of five will normally be observable ;o a user 
located anywhere on earth. 

The satellites transmit on two frequencies: LI = 1575.42 MHz and L2 = 1227.6 
MHz. The satellites transmit their signals using spread spectrum techniques employing 
two different spreading functions: a 1.023 MHz coarse/acquisition (C/A) code on LI 
only and a 10.23 MHz precision (P) code on both LI and L2. Both P-code and C/A-code 
enable a receiver to determine the range between the satellite and the user. Superim- 
posed on both the P-code and the C/A-code is the NAVIGATION message (NAV-msg), 
containing satellite ephemeris data, atmospheric propagation correction data, and satel- 
lite clock-bias information. The TANS Quadrex GPS receiver flown on STS-51 utilizes 
only the C/A-code on the LI frequency carrier. 

Two levels of navigation are provided by the GPS; these are Precise Positioning 
Service (PPS) and Standard Positioning Service (SPS). The PPS is a highly accurate 
positioning, velocity, and timing service which is made available only to authorized 
users through cryptographic keys. The SPS is a less accurate positioning and timing ser- 
vice which is available to all GPS users. The TANS Quadrex GPS receiver flown on 
STS-51 is an SPS receiver. In the future, receivers to be installed as part of the Orbiter 
GPS navigation system will be PPS units. (Kachmar, et.al., 1993, pp. 313-326) 
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The SPS is specified to provide a 100 meter (95% confidence) horizontal accuracy 
to any GPS user during peacetime. This is approximately equal to 156 meters three- 
dimensional (3-D) (95% confidence) accuracy. SPS receivers can achieve approximately 
337 nanosecond (95% confidence) Coordinated Universal Time (UTC) time transfer 
accuracy. The SPS is primarily intended for civilian purposes, although it has many 
peacetime military uses as well. The SPS horizontal accuracy specification includes the 
peacetime degradation of Selective Availability (SA) which is the dominant SPS error 
source. 1 The SA position error distribution resembles a Gaussian distribution with a 
long-term mean of zero. The SPS peacetime velocity degradation due to SA is classified. 

The ranging codes broadcast by the satellites enable a GPS receiver to measure the 
transit time of the signals and thereby determine the range between a satellite and the 
user. The NAV-msg enables a receiver to calculate the position of each satellite at the 
time of transmission of the signal. Four satellites are normally required to be simulta- 
neously “in view” of the receiver for 3-D positioning purposes. This allows the user 3-D 
position coordinates and the user clock offset to be calculated from the satellite range 
and position data. Treating the user clock offset as an unknown eliminates the require- 
ment for users to be equipped with precision clocks. Less than four satellites can be used 
if the user altitude or system time is precisely known. 

When the receiver has acquired the satellite signals from four GPS satellites, 
achieved carrier and code tracking, and has read the NAV-msg, the GPS receiver is 
ready to start navigating. The GPS receiver normally updates its pseudoranges and rela- 
tive velocities once every second. The measurements are termed pseudorange because 



1 SPS did not routinely meet accuracy specs while the GPS system was undergoing test and verifi- 
cation by the DOD, prior to December, 1993. SA effects were often varied to further degrade accuracy. 
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the clock offset of a GPS receiver introduces a bias to the true range of the satellite. The 
GPS receiver must know the GPS system time very accurately, because the satellite sig- 
nals contain the time of transmission from the satellite in GPS time. The difference in 
time between the signal leaving the satellite and arriving at the GPS receiver antenna is 
directly proportional to the range between the satellite and the GPS receiver, so it is of 
the utmost importance that the same time reference is used by both the GPS satellites and 
the GPS receiver. 

The GPS satellites carry two rubidium, and two cesium atomic frequency stan- 
dards. However, the GPS receiver is not required to have a high accuracy clock such as 
an atomic time standard. Instead, a crystal oscillator is used and the GPS receiver cor- 
rects its offset from GPS system time by making four pseudorange measurements. The 
GPS receiver can use the four pseudoranges to solve four simultaneous equations with 
four unknowns. The position equations are shown in Figure 2. When the four equations 
are solved, the GPS receiver has estimates of its position and GPS system time. The 
GPS receiver velocity is calculated using the same types of equations, using relative 
velocities instead of pseudoranges. GPS receivers perform most calculations using an 
earth-centered earth-fixed coordinate system. They then convert to an earth model 
defined by the World Geodetic System 1984 (WGS 84). WGS 84 is a very precise model 
that provides a common grid system for transformations into other coordinate systems or 
map datums. 

Satellite coverage, as measured at the user antenna, can be affected by physical 
obstructions, vehicle maneuvering or aspect, and basic receiver design. It therefore can- 
not be categorized as a GPS system requirement or specification. Coverage is defined by 
the orbits of the active satellites. The orbits determine the geometric relationships 
between the satellites and the user, which the user measures as Position Dilution of Pre- 
cision (PDOP). Since the geometric relationships continuously change as the satellites 
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Figure 2: GPS Position Equations 

move round their orbits, so does the user’s value of PDOP. PDOP is defined as the 
square root of the variances of the position errors ( PDOP = (cr x 2 + a y 2 + a z 2 )' /2 ), and 
its effects are illustrated in Figure 3. 

C. SUMMARY 

The GOOD DTO was designed to display real time GPS position and velocity data 
to the astronaut crew while on-orbit, and to record this data on hard disk for post-flight 
comparison with a ground generated Best Estimate of Trajectory (BET). It was also 
designed to display real time state vector differences between the following: 

• TANS GPS and Orbiter IMU 

• TANS GPS and SPAS GPS 

• Orbiter IMU and SPAS GPS 

The state vector differences between SPAS GPS and either the Orbiter IMU or TANS 
GPS could then be compared against on-board radar and laser rangefinder results. 
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Figure 3: Effects of PDOP 



1. State Vector Differencing 

Three prerequisites had to be satisfied before two state vectors could be com- 
pared or differenced. These were as follows: 

• Position and velocity elements were in the same frame of reference. 

• Time elements were in the same time scale. 

• Time elements were exactly matched. 

Since the Shuttle navigation system used the M50 Earth-centered Inertial 
frame of reference, and the GPS receivers utilized the WGS 84 Earth-centered Earth- 
fixed frame of reference, GPS state vectors were rotated from the WGS 84 reference 
frame to the M50 reference frame prior to comparison with the Orbiter IMU state vec- 
tor. This coordinate transformation will be addressed at length in a later chapter. 

Since the Shuttle navigation system used Greenwich Mean Time (GMT), and 
the GPS system used GPS time, GPS state vector times were adjusted to GMT prior to 
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comparison with Orbiter IMU state vectors. The Orbiter's clock was set according to the 
National Bureau of Standards UTC standard, making the GMT time scale equivalent to 
the UTC time scale in this application. The difference between UTC and GPS time was 
transmitted in the NAV-msg by GPS satellites, and this value was used for the adjust- 
ment. 

Since the state vectors being compared were produced by independent sys- 
tems, in general, the time elements did not match. The GOOD DTO utilized Cowell’s 
method to “propagate” the earlier state vector forward in time until its time element 
matched the time element of the state vector it was being compared to. Although the dif- 
ference between time elements was typically less than a second, it was significant at 
orbital velocities of several kilometers per second. Quick and accurate propagation of 
states is an active area of research, but will not be addressed further in these pages . 1 

2. State Vector Filtering 

Both the Orbiter IMU, and the SPAS GPS state vector outputs were filtered 
in order to smooth the outputs over time, and to improve accuracy. The TANS GPS stat- 
evector outputs, however, did not undergo filtering. Though use of a Kalman filter to 
smooth the output would likely improve the TANS accuracy, a filter was not imple- 
mented for the GOOD DTO due to time constraints. A Kalman filter has since been 
designed for use with the TANS and will be addressed in a later chapter. 



^ee the master’s thesis of LT Lester Makepeace for a discussion of propagation, and an alternative 
to Cowell’s method for this application. 
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II. EXPERIMENT DESCRIPTION 



A. INTRODUCTION 

The test objectives for this experiment were to demonstrate GPS on-orbit perfor- 
mance at a relatively low cost. To meet this objective, NASA needed to: 

• minimize interfaces to the Orbiter; 

• use “off-the-shelf” GPS technology; 

• designate this flight test as a non-critical (Detailed Test Objective) DTO; and, 

• limit hardware/software certification and qualification to ensuring crew safety. 

Other objectives for this flight included collecting on-orbit GPS data to be pro- 
cessed post flight, demonstrating the GPS performance to STS-51 crew real time, and 
evaluating potential future use for this hardware and software. 

The GOOD (GPS On-Orbit Demonstration) software was developed for use on a 
GRID 386 laptop computer operating at 10 MHz. The desire was to provide the crew 
with the capability to command and control the GPS receiver, and to display and record 
GPS data for real time and postflight analysis. In February 1993, when our Naval Post- 
graduate School team (LT Lee Barker, LT Les Makepeace, LT Steve Rehwald, and LT 
Carolyn Tyler) arrived at NASA, Johnson Space Center, the flight hardware had already 
been selected. The software used to interface with the TANS GPS receiver was being 
fine-tuned to NASA’s needs. Software used to provide real time analysis of the GPS data 
had not been completed. 



B. LEVEL I: HARDWARE 

The hardware selected by the NASA team consisted of the TANS receiver and its 
associated preamplifier, cabling, and three antennas. Foam spacers were used to keep 
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the antennas away from damaging the sensitive ultraviolet radiation protective coating on 
the Orbiter’s windows. Radio Frequency (RF) absorbers were placed behind the anten- 
nas to minimize the RF disturbances in the crew cabin. Finally, the entire assembly was 
attached to the window using velcro straps. As a safety requirement, a JAM (Junction 
Adapter Module) was designed and built. The JAM was the only electrical interface to 
the Orbiter and was required to protect the Orbiter from any adverse electrical behavior 
generated by the GPS experiment. The hardware set up used aboard Discovery is shown 
in Figure 4. 
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Figure 4: Hardware Setup 
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A more specific description of the TANS receiver used in this experiment is listed 



in Table 1. 

Table 1: TANS RECEIVER DESCRIPTION 



TANS Description 


Specification 


Code / Carrier Tracked 


C/A code, LI 


Channels 


6 


Antenna input signals 


up to 4 


Position Accuracy 


25 meters (SEP) without SA 




100 meters (2DRMS) with SA 


Velocity Accuracy 


0.2 m/s without SA 




classified with SA 


Time Accuracy 


1 microsecond of UTC 


Dimensions: 




- Receiver / Antenna 


5”x9.5”x2.2” / 3.75”x4”x0.54” 


Prime Power 


3.5 Watts @ 28 VDC 


Weight: Receiver / Antenna 


3.5 lbs/ 0.4 lbs 


Dynamic Capability: 




- Velocity 


8000 m/s 


- Acceleration 


4 g’. c 


- Jerk 


2 g ’s/sec 


Data Interface 


RS 422 dual channel 


Temperature: 




- Operating / Non-Operating 


-40 to 70 / -55 to 85 degrees C 


Altitude 


1 100 nautical miles 


Vibration 


0.04 g /Hz, 100 to 1100 Hz 


Shock 


40 g / 11 ms, 75 g / 6 ms 


Humidity 


100% condensing 
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The TANS receiver was chosen because it was space-configured and it met the 
main objective of being available at a relatively low cost. A disadvantage to using this 
receiver was that it had Standard Positioning Service (SPS) capability and not the pre- 
ferred Precise Positioning Service (PPS). The latter capability would have increased 
cost, slowed progress and put special restrictions on this experiment which would have 
prevented it from making STS-51’s scheduled flight deadline. In the future, NASA 
intends to use the PPS capability which improves accuracy by 100 to 156 meters as com- 
pared to the SPS receiver. Another disadvantage to using this receiver was its use of a 
deterministic point solution design instead of a filter. If the receiver had incorporated a 
Kalman Filter as a part of its design, the processor unit would perform calculations 
based on a filtering design and not based on user-selected specifications. In this flight 
test, as an example, NASA chose a 3-D Manual selection, an option the TANS provides 
the user. If the 3-D Manual is on, a three dimensional solution will not be calculated 
unless four satellites are in view and meet certain requirements. 

There were several settings made to the TANS which kept this receiver from being 
tested in its best configuration. The optimum TANS receiver/antenna performance was 
found to occur when antennas were placed on an unobstructed flat surface, looking 
straight up into the GPS constellation. Unfortunately, the Orbiter does not fly in an atti- 
tude or provide a window set up to facilitate such an optimum antenna placement, but 
rather the Shuttle flies in a left, right, or both wings down attitude (payload bay down 
towards the earth) with small restrictive windows. As one might guess, the worst posi- 
tion for the receiver was when the Orbiter was in the latter position, with all three anten- 
nas facing mostly away from the GPS satellites. Due to this major drawback, NASA 
made special exceptions to important settings like Position Dilution of Precision 
(PDOP), Signal Level Mask and Elevation Angle. 
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The PDOP mask, a user-specified setting, should normally be set to a maximum of 
six. This setting ensures that four satellites are in a satisfactory geometrical position. 
NASA set the PDOP at 20 because there was concern that a setting of six might be too 
restrictive, due to limited viewing angles and Orbiter attitudes, which would cause the 
receiver to make very few solutions. 

NASA set the Signal Level Mask specification at five AMU’s (a linearized mea- 
surement of noise, usually expressed in decibels) which is a very low setting. Nominal 
levels are in the 15 to 25 range. If this setting had been lower, more cockpit noise would 
have interfered. As it was, the PGSC’s (the other laptops in use) were picked up by the 
antennas and tracked by the receiver, even at levels below five AMU’s. Trimble, since 
this experiment, has modified the TANS software package to correct for this problem. 
The need for this low setting was apparent after looking at the way a satellite was 
tracked. In order to track the satellite as long as possible, the signal cut off had to be 
low. Also, the elevation angle was required to be low to enhance tracking, and in this 
case, NASA set the angle to zero. This encouraged error due to multipath effects and the 
tendency of the antenna to pick up competing signals (i.e. the PGSC) when tracking the 
satellites through low angles. In Chapter V, the above error inducing settings will be 
covered in greater detail. 

C. LEVEL I AND II: SOFTWARE OVERVIEW 

Software was classified as either Level I or Level II. These categories reflected not 
only their level of importance, but the minimum requirements needed for the experiment 
to fly aboard the shuttle. 

At the important but basic first level, this experiment required the ability to oper- 
ate the TANS receiver and display certain data to a PGSC (a GRID 386/10) screen. It 
had to store state vectors and engineering data to files, with the ability to downlink this 
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data. The downlink capability was important should a situation arise when the astronaut 
crew required troubleshooting assistance from ground experts during the flight. . 

At the more advanced level. Level 11, NASA desired some additional capability, 
such as the ability to input an Orbiter state vector either manually or automatically to 
compare the GPS state vector measurement to the Orbiter's. NASA also wanted to com- 
pare the TANS GPS measurements with another GPS receiver used on the STS-51 pay- 
load, ORFEUS-SPAS. One of the primary missions for STS-51 was to carry this 
German-made satellite into space, release it to operate independently for several days, 
and then rendezvous to recover it prior to returning home. The NASA engineers saw this 
as an opportunity to study the relative GPS technique. 

At a minimum, NASA wanted to collect the GPS data via the TANS receiver for 
postflight studies. Ultimately, even after both Levels I and II were fully developed, this 
GOOD test was only an experiment (or DTO). It was to be operated by the STS-51 crew 
on a not-to-interfere basis, only. An example of interference occurred during the rendez- 
vous with the ORFEUS-SPAS, which was one of the specified phases of flight to record 
information for postflight study. The TANS experiment could not be run because the 
antennas, strapped in the windows, adversely blocked the crew’s view for rendezvous 
and as a safety of flight concern interfered with the crew and their duties. 

D. LEVEL I: SPECIFICS 

Available in the Level I software were four interface displays. One display showed 
the user a current TANS configuration set-up. A second display was used to send com- 
mands and requests to the TANS receiver. A third showed the Orbiter’s location on a 
world map, and the forth display showed data for crew monitoring. This data consisted 
of six rows of information for six channels and their related channel ID, satellite ID, 
acquisition flag, ephemeris flag, azimuth, elevation, and doppler. 
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Also provided to the user was position and velocity in two different coordinate sys- 
tems. The TANS receiver made measurement calculations using WGS-84 . Its output for 
position and velocity in this frame was either in cartesian coordinates or latitude, longi- 
tude, and altitude. A non-trivial coordinate transformation was performed on the WGS- 
84 position and velocity to translate them into one of the key reference frames used by 
the Orbiter, M50. 

E. LEVEL II: SPECIFICS 

The Level II software provided real-time graphical displays of the TANS or 
ORFEUS-SPAS GPS state vector measurements and GPS/Orbiter state vector compari- 
sons. Figure 5 demonstrates the entire GPS configuration in a block diagram. 




Figure 5: Experiment Block Diagram 



In words, the TANS GPS information traveled from the receiver via the RS-422 
connection to the PGSC 715 (the laptop dedicated to GPS). The Orbiter state vector and 
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ORFEUS-SPAS GPS data traveled from the PCMMU via another GRID 1535 laptop 
running the PCDecom program through the RS-232 connection to the PGSC 715 to be 
manipulated in the Level II code. Part of the manipulation designs were to use the Orbit- 
er’s GPS (from TANS) and ORFEUS-SPAS GPS state vectors in a rendezvous program. 

Much of the Level II code was written by the NPS team, with major guidance and 
assistance from a computer programming wizard, and author of the PCDecom program, 
Mr. Tom Silva. Greater detail about the flight code and mathematical derivations is sup- 
plied in the theses written by LT Lee Barker and LT Les Makepeace. 
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m. STATE VECTOR DIFFERENCING 



A. DATA ANALYSIS 

1. Reference Trajectory 

Following the STS-51 mission, trajectory data for the GOOD DTO was 
compiled by the Shuttle Navigation group at NASA Johnson Space Center. In lieu of a 
true Best Estimate of Trajectory, a reference trajectory, generated by propagating 
between real-time navigation solutions as computed during the mission, was created. 
The resulting trajectory had an accuracy of 225 meters and 0.15 meters per second in 
total position and velocity. This accuracy was deemed sufficient given the less than opti- 
mal test conditions of the DTO. 

Orbiter state vectors, whose time elements exactly matched those of GPS 
state vectors 1 were computed by interpolating within the reference trajectory using a 
cubic spline technique. This method introduced an estimated five meters of additional 
position error into the reference data. 

2. GPS Data 

Most of the data output by the TANS GPS receiver during the GOOD DTO 
was recorded for postflight analysis. A summary of recorded data is shown in Table 2. 
Due to the sheer volume of data (over 15MB) recorded, a representative sample was 
analyzed. Of approximately 2500 recorded GPS fixes, 369 (roughly 15%) were chosen 
from three different time periods for analysis. Selected double precision XYZ, and 
Velocity XYZ ECEF binary formatted data packets were reformatted into ASCII text 
GPS state vectors using the TANSPOST program supplied by NASA. 
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TANS GPS position and velocity solution times-of-fix were always matched. 
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Table 2: GPS DATA OUTPUT 



Packet Description 


Packet Size (bytes) 


Rate 


Double Precision XYZ 


40 


0.33 Hz 


Double Precision LLA 


40 


0.33 Hz 


Velocity, XYZ ECEF 


24 


0.33 Hz 


Raw measurement output 


30-36 


0.33 Hz 


Satellite Tracking Status 


28 


0.33 Hz 


Satellite Selection 


25 


0.033 Hz 


Report Operating Parameters 


21 


once+any change 


Report Control Options 


8 


once+any change 


Health of TANS 


6 


0.033 Hz 


Almanac Information 


70 


On request 


Almanac Health Page, Time & Week 


41 


On request 


Ionospheric Parameters 


44 


On request 


UTC Parameters 


43 


On request 


Oscillator Offset 


8 


On request 


Satellite Ephemeris Status 


171 


On request 



3. Data Reduction 

Before a comparison could be made between the Orbiter reference state vec- 
tors and the TANS GPS state vectors, a coordinate rotation was required to transform 
the GPS state vectors from the WGS 84 Earth-centered Earth-fixed frame of reference to 
the M50 mean inertial frame of reference. Since the standard epoch of B1950.0, upon 
which M50 is based, was superseded in 1976 by the standard epoch of J2000.0, very lit- 
tle information appears in the current literature for working with the M50 system. 
Instead, direction is given for making transformations to and from the J2000.0 mean 
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inertial coordinate system. While NASA continues to utilize older transformation theory 
to convert Earth-centered Earth-fixed coordinates directly to M50, the authors chose to 
make the transformation first into the J2000.0 mean inertial system, and then into the 
M50 mean inertial system using a constant transformation matrix published by NASA 
Jet Propulsion Laboratory (Standish, 1982, pp. 297-302) to relate the M50 and J2000.0 
systems. The transformation into J2000.0 allowed the use of the newer 1980 IAU theory 
of nutation for improved accuracy. 

4. Coordinate Transformations 

The WGS 84 system is more accurately referred to as a geopotential model 
that has adopted the Conventional Terrestrial System (CTS) (1984.0), defined by the 
Bureau International de l’Heure, as its reference frame. Earth’s center of mass is the ori- 
gin of the CTS, as well as the M50 and J2000.0 systems. The Z-axis of the CTS is 
known as the Conventional Terrestrial Pole (CTP). The X-axis of the CTS is the Zero 
meridian, and is used to derive Universal Time, specifically UT1, in the same way the 
Greenwich meridian is used to derive GMT. The Y-axis completes a right-handed coor- 
dinate system. (DMA-TR-8350.2, 1987, p. 2-1) 

Transformation of WGS 84 coordinates into M50 coordinates requires two 
3x3 rotation matrices, each computed for the time of the state vector being transformed. 
Transformation of position vectors actually requires only a single rotation matrix. Subse- 
quent discussion will refer to this matrix as “matrix 1.” Transformation of velocity vec- 
tors, however, requires a second matrix to account for the rate of change within matrix 
1. This matrix will be referred to as “matrix 2.” 

Matrix 1 is actually the product of five separate 3x3 matrices, premultiplied 
to form a single matrix. The letters A, B, C, P, and M will be used to denote these five 
matrices. A applies two rotations for polar motion. B applies a single rotation for side- 
real time (Earth rotation). C applies three rotations for astronomic nutation. P applies 
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three rotations for general precession, and M applies three rotations for standard epoch 
conversion. All matrices with the exception of M are time varying. The M50 position 
vector is the product of the transpose of matrix 1, and the WGS 84 position vector. 

Matrix 2 is also the product of five separate 3x3 matrices. Four of the matri- 
ces used to compute matrix l. A, C, P, and M, are also used for matrix 2. Unlike the 
other time varying matrices, the rate of change of the B matrix is significant; so another 
matrix, denoted by B, replaces the B matrix in matrix 2. B is the rate of change of the B 
matrix. The M50 velocity vector is computed in three steps. In step one, the product of 
the transpose of matrix 2 and the WGS 84 position vector is computed. In step two, the 
product of the transpose of matrix 1 and the WGS 84 velocity vector is computed. In 
step three, the two resulting vectors are added vectorially to form the M50 velocity vec- 
tor. The transformation methodology for position coordinates is shown in Equation 1 
and the transformation methodology for velocity coordinates is shown in Equation 2. 




1 

= (aBCPm] [xyz] 



T 

WGS-84 




T 

WGS-84 



(Eq: 1) 



(Eq; 2) 



Methods for computing the A, B, B , C, and P matrices were taken from The 
Astronomical Almanac, The Explanatory Supplement to the Astronomical Almanac, and 
Defense Mapping Agency Technical Report 8350.2. These methods are summarized in 
sub-paragraphs a through d. 

a. Polar Motion 

Polar motion parameters, x p and y p , for the dates encompassing the 
STS-51 mission, were obtained from the U.S. Naval Observatory. The A matrix consists 
of a rotation about the Y-axis by angle -r p and a rotation about the X-axis by angle -y p . 
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The maximum amplitude of these parameters is approximately 0.3 arc seconds. This 
corresponds to about 10 meters of position difference at Shuttle altitudes. 
b. Sidereal Time 

The B matrix consists of one rotation about the Z-axis by an angle of 
A, where A = H 0 + AH +co*( t - At ). H 0 is Greenwich Mean Sidereal Time (GMST) at 
0 h UT1 on the day of interest. Since 1984, GMST has been defined by Equation 3, 
where T u is the number of centuries elapsed since 12 h UT1 on 2000 January 1. The 
result is in units of seconds of sidereal time, and may be converted to arc on the basis of 
one revolution per 24 hours of sidereal time. 

GMST1 of0 h UTl = 24110.54841 + 8640 184.8 12866 T + 0.0931047^ -6.2x 10" 6 7^ 

(Eq:3) 

AH is the Equation of the Equinoxes. It equals arctan( cose tanAy ), 
where e is true obliquity of the ecliptic, and A \\r is nutation in longitude. Both e and A\|/ 
are computed in the course of generating the C matrix, to* is the Earth rotation rate in a 
precessing reference frame, and is equal to (o'-t- m, where (o' is the Earth’s inertial rota- 
tion rate, a constant, and m is equal to 7.086 x 10‘ 12 + 4.3 x 10‘ 15 T U . Time t is the time 
of the state vector being transformed in seconds since the beginning of the day UTC, and 
At is the difference between UT1 and UTC. Values of UT1 minus UTC were obtained 
from the U.S. Naval Observatory for the dates encompassing the STS-51 mission. These 
values are kept below 0.7 seconds through introduction of leap seconds into UTC. One 
second of time corresponds to 15 arc seconds, or about 485 meters of position difference 
at Shuttle altitudes. 

(1) Change in Sidereal Time Matrix. The B matrix is defined as 
shown in Equation 4. 
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